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The scattering of coherent X-rays from dynamically evolving systems is 
currently becoming experimentally feasible. The scattered beam produces a 
pattern of bright and dark speckles, which fluctuate almost independently in 
time and can be used to study the dynamics of the system. Here we report 
large-scale computer simulations of the speckle dynamics for a phase ordering 
system, using a two-dimensional model quenched through an order-disorder 
transition into the two-phase regime. The intensity at each wave vector k 
is found to be an exponentially distributed random variable. The scaling 
hypothesis is extended to the two-time correlation function of the scattering 
intensity at a given wave vector, Corr(k; t\, £2)- The characteristic decay time 
difference for the correlation function, \t±— t2\ c , is found to scale as k\ti+t2\ 1 ^ 2 ■ 
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I. INTRODUCTION 



When electromagnetic radiation is scattered from an inhomogeneous material, it un- 
dergoes phase changes which depend on which part of the material it scatters from. The 
resulting pattern of constructive and destructive interference at the detector depends sensi- 
tively on the geometry of domains in the scattering region. If the incident radiation is not 
coherent over the scattering region, this dependence is washed out by the phase differences 
in the incident beam. However, if the incident beam can be made coherent over the scat- 
tering volume, a speckled scattering pattern such as that shown in Fig. 1 is observed. The 
intensity fluctuations associated with the speckle pattern are the basis of photon correlation 
experiments. Over the last three decades, lasers have become available as sources of coherent 
radiation at wavelengths ranging from the ultraviolet to the infrared [0. 

The much shorter wavelength of X-rays allows materials to be studied at smaller length 
scales than is possible with light, and the greater penetration allows the study of optically 
opaque materials. However, coherent X-rays have only become readily available with the 
emergence of high-intensity synchrotron photon sources. Recently, speckle patterns have 
been observed in several experiments using coherent synchrotron X-rays. For example, 
diffusion rates in gold colloids were determined from the time Brownian motion requires to 
change the X-ray speckle pattern [Q,^J , and speckle from coherent X-rays has also been used 
to study equilibrium fluctuations in FesAl near an order-disorder transition Q. 

An important situation for analyzing phase ordering processes is when the system is 
characterized by a single length scale R, corresponding to the average domain size, that 
grows asymptotically with time t as R cx t n , where n is known as the growth exponent. 
If the dynamics are controlled by diffusion, n — 1/2 for a non-conserved order parameter 
and n = 1/3 for a conserved order parameter. The significance of this for scattering experi- 
ments is that the scattering intensity at wavevector k is proportional to the structure factor 
S(k,t). At all times in the asymptotic scaling regime the average structure factor, (S(k,t)), 
can be expressed in terms of a scaling function which depends on k and t only through a 
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dimensionless wavevector q oc kt n \ 

t~ dn (S(k,t))cx(I(q)} , (1) 

where d is the spatial dimensionality. 

An alternative way to consider this scaling behavior, which is more useful in describing 
the time dependence of the fluctuations in the scattering intensities at different wavevectors, 
is the following. One observes that the dynamics of structures on different length scales 
k~ l become statistically indistinguishable when viewed on the corresponding time scales 
9 oc tk x l n = q x l n . This can be expressed in a scaling relation equivalent to Eq. ([!]): 

k d (S{k,t)) oc (F(6)) = 9 nd (I{6 n )) . (2) 

In this paper we present a numerical study of a simple model for the scattering speckle 
observed in coherent X-ray scattering from systems undergoing phase ordering following 
a quench through a second-order order-disorder phase transition. This is the standard 
time-dependent Ginzburg-Landau (TDGL) equation with non-conserved order parameter, 
equivalent to Model A in the classification scheme of Hohenberg and Halperin ||. Using 
standard Langevin dynamics, we study the fluctuations in k-space and time of the structure 
factor. A representative scattering pattern for the simulations conducted here is presented 
in Fig. 1 as a log-scale density plot. Because the simulation samples the order parameter at 
L d locations in real space, the scattering function can only be sampled at L d independent 
locations in k-space. This corresponds in a natural way to the independent speckles seen in 
experiments. 

Based on the scaling relations in Eqs. flU) and fl2|) one might naively expect that the 
position of an individual speckle should move as k oc t~ n as the average structure factor 
changes according to the scaling form. This is not the case. Instead, individual speckles 
are observed to brighten and dim with significant fluctuations around the time-dependent 
average intensity. The evolution of the simulated speckle intensities for a cut along the k x 
axis for a single simulation is presented in Fig. 2. They bear a striking similarity to recent 
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experimental data for CU3A11 |J. From this figure it is clear that as the system evolves, 
individual speckles do not travel through k-space. In the simulations reported here, the 
correlations between neighboring speckles were found to be on the order of 1%, so each 
speckle could be treated as independent. 

In this paper the time-time correlation in the scattering intensity at fixed k is obtained 
numerically via the time-dependent Ginzburg-Landau model for a non-conserved order pa- 
rameter. The details of this approach are presented in the next section. Section III presents 
the results. The scattering intensity is shown to be an exponentially distributed random 
variable. The simulation results are also used to test the scaling ansatz for the two-time 
intensity correlation. The results here suggest that the intensity correlation time grows as 
kt 1 / 2 . Section IV presents a review of this work along with some future directions we intend 
to pursue. 



II. METHOD 

We performed a standard simulation of the dynamics of phase ordering following a quench 
through an order-disorder phase transition in a system described by a nonconserved scalar 
order-parameter field ^(x, t). We used the time-dependent Ginzburg-Landau equation, 

#(x,t) W 
which we here have expressed in dimensionless variables for convenience. The first term in 
Eq. (H) corresponds to deterministic relaxation towards a minimum value of the free energy 
jF[^(x, £)], and the second term represents thermal noise. Since the late-time dynamics 
are controlled by a zero-temperature fixed point, we ignore the noise term, using only the 
deterministic part of Eq. (0) • 

We employed the standard Ginzburg-Landau- Wilson ip 4 free energy || with unit param- 
eters. The deterministic dynamical equation then becomes 

^M=(l + V>(x,t)-^ 3 (x,t). (4) 
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Note that by choosing the customary parameters in the free energy and the dynamical equa- 
tion equal to unity, we have merely specified the length, concetration, and time scales in the 
simulation. Equation (^) can therefore be used to represent the low-temperature dynamics 
of Model A without loss of generality. Since thermal fluctuations are explicitly excluded, 
the randomness they introduce into the domain pattern at early times J7J is implemented 
by an initial condition where ^(x, 0) consists of independent Gaussian random numbers 
with mean and standard deviation 0.1. The simulations were conducted on square lat- 
tices with periodic boundary conditions, lattice constant Ax = 1 , and a system size of 
L x = L y = L = 1024. The Laplacian in Eq. (f|) was approximated by the standard four- 
neighbor discretization, and a simple Euler integration scheme with At = 0.05 was used to 
collect data every 20 time units up to a maximum of t = 1000. Results were averaged over 
10 independent initial conditions. 

The Fourier transform ip of ip is defined as 

^(k,t) = -L=^^(x,t)e ik - x , (5) 

Vij X 



where the fact that the lattice spacing is unity in all directions has been used. The Brillouin 
zone is defined by the discrete set of wavevectors, k x , k y = 2nj/L with j £ {0, ±1, ±2, . . . ± 
(L/2 — 1), L/2}. The structure factor of the system is 

5(k,t) = |^(k,t)| 2 (6) 

and is proportional to the scattering intensity observed experimentally. To be consistent 
with the numerical integration, the magnitude of the wavevector, k(k), is defined using the 
operator relation 

- [k(k)} 2 ^(k, t) = -L£ e* k - x V 2 ^(x, t) • (7) 
V L a x 



Substituting the discrete version of the Laplacian, one gets 

1/2 



fc(k) 



2 I d — cos (k c 



a=l 



(8) 
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This has been used for the two-time correlation function discussed in the next section where 
the value k is used to bin data during circular averaging and to determine the value of 
scaling variables. Because of lattice effects we consider only those wavevectors with k < 1, 
for which the finite interface thickness can be ignored in this work. 

III. RESULTS 

The measured intensity at a given (k, t) is a random variable with characteristics deter- 
mined by the scattering system. As long as the independent domains are randomly located 
and oriented, and the domain size is much smaller than the system size, one may expect the 
central limit theorem to imply that the Fourier transformed order parameter, ^(k, t), is well 
approximated by a complex Gaussian random field, even if the real-space order-parameter 
field, ip(x,t), is not Gaussian. One consequence of this is that the probability density func- 
tion for the scattering intensity S(h, t) can be shown to be exponential for all k and t 
The probability density for the normalized intensities, s(k, t) = S(k, t) / (S(k, £)), is therefore 

P(s) = exp (-s) , (9) 

independent of k and t. The probability density P(s) is normalized and has unit mean 
and standard deviation. Since P(s) is identical for all values of k and t, only one density 
function needs to be constructed. This leads to very good statistics because on the order 
of 10 7 samples are available. The results for all k and t, with (S(k,t)) found by averaging 
over wavevectors with equivalent values of k 2 x + ky, are presented in log-linear form in Fig. 3. 
The solid line is the expected density P(s) = exp(— s). For the simulation the probability 
density is exponential for all x ^ 12, although deviations may occur above this. We have 
verified that the distribution at given values of k and t also is exponential. 

One important consequence of an exponentially distributed intensity is that the average 
scattered intensity at (k, t) is equal to the standard deviation of that intensity. This property 
can be used to simplify the normalized correlation function 
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Corr (Ml ,t 2 ) = ,8^)8^)) -(S^jS^ 



5(k,t 1 )5(k,t 2 ))-(5(k,t 1 ))(5(k,t 2 ; 



's(k, tl )){s(k,t 2 )) ' 

which has a value between —1 and 1 by construction. At t\ = t 2 , the same-time limit, 
the exponential nature of s makes the correlation unity as one would expect. As the two 
measurement times become widely separated, the values of the intensity become independent 
and the correlation decays to zero. For this relaxational system negative correlations are not 
expected. The scaling ansatz extended to this situation allows comparison of the scattering 
at different k and t by 



F(9 1 )F(9 2 ) 

Corr(0i, 9 2 ) = / . . ' - 1 , (12) 



where #i j2 = 2/c 2 t 12 are the rescaled times from Eq. (fj). The ansatz has been tested with 
the simulation data. Because of the amount of data involved, it was necessary to bin the 
results. The results for (F(9)) were binned by t and 9. For each pair of times, the value 
of (F(9i)F(9 2 )) for a single k was found and then placed into a bin in the (9\,9 2 ) plane. 
The average result for each bin was later found. The loss of precision from this procedure 
should not significantly affect the results. This function is most easily viewed as a contour 
plot in the (9i, 9 2 ) plane, as in Fig. 4. In this figure a Corr($i, 9 2 ) = 1 contour, which is not 
shown, extends along the diagonal. Moving away from that line, contours at values of 0.5, 
0.1, and 0.025 are shown. Scatter in the data becomes dominant for values less than this. 
The striking feature of this figure is that the correlation function becomes "broader" for 
later times, i.e. the correlation time for the speckle intensity becomes longer as the phase 
ordering proceeds. 

An alternative set of variables for this problem are r = (9\ + 2 )/2 and 5 = 9 2 - 9±. 
A constant value of r corresponds to a line perpendicular to the 9\ = 9 2 diagonal, while 5 
measures the distance (in units of scaled time) away from the diagonal. For a given value 
of r, one can measure the characteristic time difference, 8 C , required for the correlation to 



decay to a value of 1/2. The values of S c estimated by linear interpolation from the simulation 
results are presented in log- log form in Fig. 5. For large values of r, the relationship appears 
to be a power law. The solid line represents a least-squares fit which gives an exponent 
of 0.507 ± 0.002. Analytic arguments indicating that the exponent is exactly 1/2 will be 
presented elsewhere ||. Furthermore, at these r the correlation as a function of 5 can be 
collapsed onto a single master curve using S c . The results for a wide range of r are presented 
in Fig. 6, where the interpolated value of S c has been used. The collapse of the data is quite 
good, indicating that at large r the correlations can indeed be described by a single curve. 



IV. CONCLUSIONS 

Intensity fluctuations in the scattering of coherent light from a two-dimensional system 
undergoing an order- disorder transition have been investigated by numerical simulation. The 
ratio of the measured intensity at given (k, t) to its mean value at the same (k, t) is found to 
be an exponentially distributed random number for the dynamically evolving system. The 
two-time intensity correlation function is found to obey scaling. At late times, its shape in 
the direction perpendicular to the t\ = t 2 line, in units of k 2 (t 2 — ti), is constant except for 
a characteristic time difference which increases as the system evolves. This characteristic 
decay time is found to scale as k(t\ + 1 2 ) 1 ^ 2 . 

Recently we have found that when the domain size is much smaller than the system 
size, the speckle correlation function defined in Eq. ( |10"D can be related to the two-time 
correlation function of the order parameter in real space. Theories for this latter function 



exist [ 10 , 1 1 1 . We are currently analyzing these theories in terms of the predictions made 



above and investigating how coherent X-ray scattering could be used to test these theories 
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FIGURE CAPTIONS 



Fig. 1 Example of the speckled scattering intensity from one measurement of the simulation. 
This is a log-scale intensity plot, with lighter shades indicating brighter speckles. The 
k = origin is at the center of the figure. Speckles do not shift in k-space, but their 
intensities fluctuate strongly around the (k, t) dependent average value. 

Fig. 2 Example of the evolution of certain speckles with time. This is a cut along the k y = 
axis for one simulation; the intensity scale matches that for Fig. 1. The earliest time 
is at the bottom of the figure and time increases to the latest simulation time at the 
top. k x = is the left hand column and k x increases towards the right. From this 
figure it is clear that individual speckles do not move in k-space and that neighboring 
speckles are relatively uncorrelated. The fluctuations in the speckle intensity are more 
apparent at higher wavevectors. 

Fig. 3 The probability density for the normalized speckle intensity, s = S(k,t)/(S(k,t)). 
The solid line is the theoretical density of an exponentially distributed random variable. 

Fig. 4 Contour plot of the scaled two-time correlation function. The correlation at 9\ = 9 2 is 
unity by construction and is not shown. The contours moving away from this diagonal 
are at 0.5, 0.1, and 0.25. The figure shows that the speckle intensity stays correlated 
longer as Q\ + 9 2 increases. 

Fig. 5 The characteristic decay time difference, 5 C , as a function of r oc t\ + t 2 . For larger 
values of r, the relationship is seen to be a power law. The solid line is a least-squares 
fit which finds 5 C ~ r °- 507±a002 . 

Fig. 6 Collapse of the intensity-intensity correlation for different values of r oc t\ + t 2 , where 
5 oc t 2 — ti and 5 C is shown in Fig. 5. The figure shows that the correlation can be 
collapsed onto a single curve for a large range of r. 
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